Two-Stage Formation Model and Helicity of Gold Nanowires 
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A model for formation of helical multishell gold nanowires is proposed and is confirmed with the 
quantum mechanical molecular dynamics simulations. The model can explain the magic number of 
the helical gold nanowires in the multishell structure. The reconstruction from ideal non-helical to 
realistic helical nanowires consists of two stages: dissociations of atoms on the outermost shell from 
atoms on the inner shell and slip deformations of atom rows generating (lll)-like structure on the 
outermost shell. The elementary processes are governed by competition between energy loss and 
gain by s- and d-electrons together with the width of the d-band. The possibility for the helical 
nanowires of platinum, silver and copper is discussed. 
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Gold nanowires synthesized by electron beam tech- 
nique have helical multishell structures. [H, 0] They are 
formed along the original [110] axis. The outermost shell 
is a (lll)-like atomic sheet and [110] atomic rows are 
helical around the nanowire axis. The histogram of ap- 
pearance of the diameters of nanowires has peaks of the 
numbers of atoms discretely at seven, eleven, thirteen, 
fourteen and so on. The difference of numbers of atoms 
between the outermost and the next outermost shells is 
seven, called "magic number", except cases of five and 
seven atoms on the outermost shell. Platinum nanowires 
were also reported to have the same type of helicity. Q 
Therefore, the helical multishell structure of nanowires 
might be generic not only for gold but also for a certain 
class of metals. 

A structure of nanowires should be determined gener- 
ally by competition between bulk and surface regions. A 
pioneering theoretical work ^ predicted, with classical 
potentials of Al and Pb, that such nanowires can form 
various noncrystalline 'exotic' structures including heli- 
cal ones. There are several theoretical works on [110] 
helical multishell gold nanowires, based on quantum me- 
chanical simulation. 0, 0, 0] Tossatti et al. reported, 
with the first principles calculations, that the tension of 
nanowires gives the minimum values when the number of 
atoms of the lateral atomic row on the outermost shell is 
seven and the nanowire is helical. Q They also showed 
that the tension does not have the minimum in model 
silver nanowires. Senger et al. achieved the first prin- 
ciples electronic structure calculations of nanowires with 
atoms from three to five on the lateral atomic row on the 
outermost shell, and showed that helical nanowires are 
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not the configuration of the minimum energy but of the 
minimum tension. Q 

In spite of these intensive studies, no theories so far ex- 
plain the formation process of the helicity and the magic 
number of helical nanowires of gold and platinum. In 
the present Letter, wc propose a two-stage model of for- 
mation of multishell helical nanowires and explain the 
origins of the helicity and the magic number. 

In the ideal [110] fee nanowire of stacking (110) sec- 
tions, there are atom rows parallel to the [110] axis and 
we can find multishell structure with no helicity. Fig- 
ure [TJ^ a) depicts three structures of the ideal fee [110] 
nanowires, 6-1, 10-4, and 12-6, where the index stands 
for the number of atoms on the lateral atom row on the 
outermost shell, that on the next outermost shell and so 
on. The structure of these ideal model nanowires satisfies 
two conditions: (a) there is no acute angle on the surface 
because of diminishing surface tension, and (b) there is 
no (001) side longer than any (111) side since the surface 
energy of a (001) surface is higher than that of (111). 

In the present model, the driving force for the helicity 
is the atom row slip. The helical structure is introduced 
when the outermost shell transforms into a folded (111)- 
like sheet. The transformation consists of two stages. 
Firstly the outermost shell should be dissociated from 
the inner shell to move freely and atom rows along [110] 
might be inserted into the outermost shell (Stage 1). Sec- 
ondly, since a (111) sheet is energetically favorable on the 
outermost shell, an atom row slips and (001) faces trans- 
form into (111) ones (Stage 2). 

The outermost shell has six more "bonds" on the lat- 
eral row than the inner shell. See 6-1, 10-4 and 12-6-1 
nanowires, shown by bold lines in Fig.[TJa). We add one 
another atom row on the outermost shell as in Fig.[TJb), 
then the outermost shell has seven more atoms on a lat- 
eral atom row than the inner shell and the outermost 
shell can have room for atom row slip. This is the origin 
of the "magic number" . When the number of atoms on 
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FIG. 2: The sections of gold nanowires: (a) 7-1 at the initial 
(left), 400 (middle) and 5,000 MD steps (right), (b) 11-4 at 
the initial (left), 400 (middle) and 6,000 MD steps (right), 
(c) 12-6-1 at the initial (left), 7,500 (middle) and 13,000 MD 
steps (right), (d) 13-6-1, at the initial (left), 750 (middle) and 
9,000 MD steps (right), (e) 14-7-1, at the initial (left), 400 
(middle) and 7,000 MD steps (right), and (f) 15-8-1 at the 
initial (left), 400 (middle) and 5,000 MD steps (right). An 
atom row slips on the outermost shell and helicity is intro- 
duced. Temperatures in the simulations are 600 K in (a), (b), 
(e), and (f) and 900 K in (c) and (d). 



FIG. 1: (a) The ideal [110] fee nanowire is intrinsically a non- 
helical multishell nanowire. (b) Expanded a surface of the 
7-1 multishell nanowire before (upper right) and after (lower 
right) surface reconstruction. Dashed lines connect the same 
atoms at the right and the left ends. 



the lateral row in the outermost shell is odd, the surface 
reconstruction brings the helicity to nanowires inevitably. 
Figure[ljb) shows the case of forming 7-1 helical nanowire 
where a slip of atom row transforms the surface region 
B, initially a (lOO)-like surface, into a (lll)-like surface 
and a folded (lll)-like surface brings the helicity. Helical 
nanowires of 7-1, 11-4, 14-7-1 and 15-8-1 are experimen- 
tally observed. [l| 

We performed, to verify the above model, the MD sim- 
ulations of gold nanowires with the tight-binging (TB) 
Hamiltonian. These TB parameters are prepared 

to represent electronic structures of bulk solids, surfaces, 
stacking faults and point defects, and have been tested in 



several other simulations for gold nanowires. 
We stack ideal (110) sections of the fee lattice, as an ini- 
tial state, as depicted in Fig. [21 nine layers in the cases 
of 7-1 and 11-4, seven layers in the cases of 12-6-1 and 
13-6-1, and six layers in the cases of 14-7-1, and 15-8-1. 
Then the numbers of atoms are from 76 to 156. We adopt 
the boundary condition of fixing the center of gravity of 
both end layers of the nanowires, and impose no external 
force. One MD step is 1 fs. Figure [2] shows the time evo- 
lution of typical sections in gold nanowires. The outer- 
most shell becomes a smooth surface and some atoms on 
it dissociate from atoms on the inner shell. The surface 



reconstruction is always observed, and no significant dif- 
ference is observed at two different temperatures, 600 K 
and 900 K. In the case of 12-6-1 nanowire (Fig. (H^c)), no 
surface reconstruction is observed at 600 K. 

A typical time evolution of MD simulations is shown 
in Fig. [3] for an 11-4 nanowire (Tig. [^b)) at 600 K. See 
also the supplementary movie. [l3| The initial state is a 
10-4 nanowire with an additional atom row, and the to- 
tal number of atoms is 143. The total energy decreases 
almost monotonically after 1,000 MD steps. The atom A 
in Fig. [31 which is on the (111) sheet and also its nearest 
neighbor atoms (e.g. B) on the (100) sheet, dissociates 
from the inner sheU before 500 MD step (Stage 1). ^ 
Then the surface atoms can move independently from 
inner sheU. From 2,000 to 5,000 MD steps, (001) sheet 
reconstructs into hexagonal (lll)-like surface with a slip 
deformation, and the helical structure appears (Stage 2). 
The inner shell rotates at the same time. These behav- 
iors with two stages agree with the proposed model. It 
should be noticed that the time evolution in the MD sim- 
ulation does not totally correspond to a real formation 
process of gold nanowires. But we believe that the real 
formation process has two different stages of short time 
process (Stage 1) and long time process (Stage 2). 

We are focusing on the change of the electronic struc- 
tures to know why the atom A dissociates from the inner 
shell in Stage 1, and why the slip deformation appears 
around the atom B in Stage 2. The local electronic den- 
sity of states (LDOS) is shown in the local coordinate 
system in which the x-axis is the nanowire axis, [110], 
the y-axis is along the lateral direction and the z-axis 
is perpendicular to the surface, so that the xy-orbital 
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FIG. 3: Change of the total energy and the atomic configu- 
ration during the formation process of the 11-4 nanowire at 
600 K. Red, blue, purple, and green atoms are named A, B, 
C, and D respectively. Within 500 MD steps, the atom A was 
disassociated from inner shell. Between 2,000 and 5,000 MD 
steps, the slip between the atom B and D brought the surface 
reconstruction, and (001) face becomes (111) face. 

should be on the surface on the nanowire. 

First we discuss the change of LDOS of the atoms A 
and C during the dissociation of the atom A (Stage 1). 

The dissociation of the atom A from the inner atoms 
changes the electronic structure of the atom A in Stage 
1 as follows. The local energy of s-orbital increases be- 
cause of reduction of its coordination number (Fig. [H^a) 
top). The local energy in yz-orbital of the atom A (not 
shown) increases since the orbital extents perpendicu- 
larly to the crystal surface and the nearest neighbor dis- 
tance increases. The local energy in xy-orbital decreases, 
since the orbital can expand more to another (111) sheet 
through the atom C because of flattening two (111) sheets 
and the band width becomes wider (Fig. IHa) middle). 
Then the energy loss and gain of the d-orbitals almost 
cancel with each other. The net loss of the local energy 
of the atom A is mainly due to energy loss of the s-orbital 
(Fig. Ufa) bottom). 

Figure m (b) shows LDOS of the atom C, the nearest 
neighbor of the atom A. There is relatively small energy 
loss of yz-orbital in Stage 1 since the coordination num- 
ber of the atom C does not change. The local energy of 
xy-orbital decreases (Fig. SJ^b) middle). Therefore, the 
local energy of the d-orbitals decreases on the atom C. 
This energy gain of d-orbitals can be attributed to the 
flattened surface structure around the atom C after the 
dissociation between the atom A and the inner shell. The 
energy of s-orbital of the atom C also decreases apprecia- 
bly (Fig. Hl^b) top) but it may be not associated primarily 
with the dissociation since the s-orbital does not always 
favor the flatter atomic configurations. 

The dissociation between the atom A and the inner 
shell in Stage 1 is governed by competition between the 



s-orbital energy loss of the atom A with reducing the co- 
ordination number and the energy gain of the d-orbitals 
due to the flattened surface. The change in the energy 
of d-orbitals is almost proportional to the d-band width, 
and gold has wider d-band width than that of copper 
and silver. Therefore, the bond dissociation of Stage 1 
is preferable for gold but not so for silver and copper. 
In fact, the MD simulation in copper does not show the 
bond dissociations of Stage 1 within 500 MD steps. 

The electronic structure of the atom B changes as fol- 
lows during the slip deformation (Stage 2). LDOS of the 
atom B is depicted in Fig. [4] (c). The local energy of the 
s-orbital decreases since the coordination number of the 
atom B increases from seven to eight (Fig. 2] (c) top). 
The local energy in the xy-orbital of the atom B also de- 
creases since the surface transforms to (lll)-like and is 
flattened (Fig.[4](c) middle). Here, the slip deformation 
is essential since it widens the area of the (lll)-like trian- 
gular surface and LDOS of xy orbital transfers its weight 
from the anti-bonding region (the high energy region) 
to the bonding region (the low energy region) (Figs. [4] 
(a), (b) and (c) middle). Therefore, both orbitals s and 
xy of the atom B result in decrease of the local energy. 
The accomplishment of slip deformation is, however, de- 
pends on the details of surrounding atoms. The atom A 
connecting with the atom B can move relatively freely 
from the inner shell since the bondings with the inner 
shell are dissociated after Stage 1. Therefore, the atom 
B slips easily and creates the helicity since the atom B 
can trail other atoms without dissociating their bonds. 
Higher temperatures were required for reconstruction of 
(001) surface in copper nanowires since the atom A in 
copper does not dissociate from the inner shell. 

The two-stage formation model holds also in 12-6-1 
and 13-6-1 nanowires. Figures [2{c) and (d) show how 
inserted atom row would behave. See supplementary 
movie, [l^ Both nanowires show surface reconstruction, 
and the number of atoms in the outermost shell is 13. 
Therefore, the inserted atoms are supplied possibly from 
outer and inner shells and positions of supplied atoms is 
not seriously crucial to the surface reconstruction, or the 
atom row slip. 

In all cases, the bond dissociation between the outer- 
most atoms and the inner shell starts near a cross edge 
of two ideal faces, since the energy gain is maximized by 
flattening a surface on the edges. The ratio of the en- 
ergy gain to the total energy due to the dissociation be- 
comes smaller in thicker nanowires and then the thicker 
nanowires do not transform. 

MD simulation for 11-4 and 12-6-1 nanowires in copper 
does not show the bond dissociation from the inner shell 
of Stage 1 at 600 K within 500 MD steps, since copper 
has narrower d-band width. At 900 K the atom row slips 
on the (001) sheet in copper nanowires, since the atom 
have dissociated from the inner shell at that temperature. 

The present analysis shows that the mechanism in both 
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FIG. 4: Partial densities of states and local energies in or- 
thogonalized basis. Partial densities of states of (a) the atom 
A, (b) the atom C, and (c) the atom B. Red and blue lines 
are at the initial state and at 500MD steps respectively in 
(a) and (b), and at 500MD steps and 5,000MD steps in (c). 
Upper, middle, and lower figures are the partial densities of 
states of s-orbital, xy-orhital, and the total density of states 
respectively. The energy changes of the densities of states are 
written at the upper right of each figure. Local coordinate 
system for each atom is written in text. Notice that, in all 
cases, the asymmetry of LDOS of xy orbital is enhanced. 



Stage 1 and Stage 2 is governed by the d-band width and 
the helical nanowires appear only among metals with a 
wider d-band. The d-band width in platinum and gold 
is commonly wider than that in lighter elements, Ag 
and Cu, and the present theory explains why platinum 
nanowire can be also formed with helicity. 

Finally, we point out that the present mechanism is 
inherent not only with nanowires but also with the bulk 
surface, e.g. the Au (001) surface reconstruction to Au 
(OOl)-hex, hexagonal close-packed one. 15, This also 
implies a reason why (001) surface reconstruction is ob- 
served in metals with 5d orbitals such as Au, Pt, and Ir, 
but not observed in metals with the outer 3d or 4d shell 
such as Ag and Cu. 

In summary, we have proposed the two-stage model 
of formation of gold multishell helical nanowires and its 
'magic number'. The model is confirmed by molecu- 
lar dynamics simulation with electronic structure and its 
mechanism is explained as a nanometer-scale size effect 



dominated by d-orbitals on nanowire surface. Helicity is 
introduced by the surface reconstruction or the slip of 
an atom row on the (001) sheet, because the triangular 
(lll)-like sheet is more preferable for d-orbitals extend- 
ing over the surface. 
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